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Abstract. A semilinear relation is a finite union of finite intersections of open and 
closed half-spaces over, for instance, the reals, the rationals, or the integers. Semilinear 
relations have been studied in connection with algebraic geometry, automata theory, and 
spatiotemporal reasoning. We consider semilinear relations over the rationals and the 
reals. Under this assumption, the computational complexity of the constraint satisfaction 
problem (CSP) is known for all finite sets containing R+ = {{x,y,z) \ x + y = z}, <, 
and {!}. These problems correspond to expansions of the linear programming feasibility 
problem. We generalise this result and fully determine the complexity for all finite sets 
of semilinear relations containing R+. This is accomplished in part by introducing an 
algorithm, based on computing affine hulls, which solves a new class of semilinear CSPs 
in polynomial time. We further analyse the complexity of linear optimisation over the 
solution set and the existence of integer solutions. 
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1. Introduction 

We work over a ground set (or domain) X, which will either be the rationals, Q, or 
the reals, M. We say that a relation R C is semilinear if it can be represented as a 
finite union of finite intersections of open and closed half-spaces in X^. Alternatively, R is 
semilinear if it is first-order definable in {R+, <, {!}} where = {(x, y, z) G X^ \ x + y = 
z} m- Semilinear relations appear in many different contexts within mathematics and 
computer science: they are, for instance, frequently encountered in algebraic geometry, 
automata theory, spatiotemporal reasoning, and computer algebra. Semilinear relations 
have also attained a fair amount of attention in connection with constraint satisfaction 
problems (CSPs). In a CSP, we are given a set of variables, a (finite or infinite) domain 
of values, and a finite set of constraints. The question is whether or not we can assign 
values to the variables so that all constraints are satisfied. From a complexity theoretical 
viewpoint, solving general constraint satisfaction problems is obviously a hard problem. 
Various ways of refining the problem can be adopted to allow a more meaningful analysis. 
A common refinement is that of introducing a constraint language] a finite set T of allowed 
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relations. One then considers the problem CSP(r) in which all constraint in the input must 
be members of F. This parameterisation of constraint satisfaction problems has proved 
to be very fruitful for CSPs over both finite and infinite domains. Since P is finite, the 
computational complexity of such a problem does not depend on the actual representation 
of the constraints. 

The complexity of finite-domain CSPs has been studied for a long time and a powerful 
algebraic toolkit has gradually formed [7| . Much of this work has been devoted to the Feder- 
Vardi conjecture naiii] which posits that every finite-domain CSP is either polynomial¬ 
time solvable or NP-complete. Infinite-domain CSPs, on the other hand, constitute a 
much more diverse set of problems. In fact, every computational problem is polynomial¬ 
time equivalent to an infinite-domain CSP [T]. Obtaining a full understanding of their 
computational complexity is thus out of the question, and some further restriction is 
necessary. In this article, this restriction will be to study semilinear relations and constraint 
languages. 

A relation R C is said to be essentially eonvex if for all p,q ^ R there are only 
finitely many points on the line segment between p and q that are not in R. A constraint 
language P is said to be essentially convex if every relation in P is essentially convex. The 
main motivation for this study is the following result; 

Theorem 1 (Bodirsky et al. [2]). Let F be a finite set of semilinear relations over Q or 
M such that {i?+, <, {!}} C P. Then, 

(1) CSP{T) is polynomial-time solvable if T is essentially eonvex, and NP-eomplete 
otherwise; and 

(2) the problem of optimizing a linear polynomial over the solution set of CSPfT) is 
polynomial-time solvable if and only if CSP(T) is polynomial-time solvable (and 
NP-hard otherwise). 

One may suspect that there are semilinear constraint languages P such that CSP(r) € P 
but P is not essentially convex. This is indeed true and we identify two such cases. In the 
first case, we consider relations with large “cavities”. It is not surprising that the algorithm 
for essentially convex relations (and the ideas behind it) cannot be applied in the presence 
of such highly non-convex relations. Thus, we introduce a new algorithm which solves 
CSPs of this type in polynomial time. It is based on computing affine hulls and the idea 
of improving an easily representable upper bound on the solution space by looking at one 
constraint at a time; a form of “local consistency” method. In the second case, we consider 
relations R that are not necessarily essentially convex, but look essentially convex when 
viewed form the origin. That is, any points p and q that witnesses a not essentially convex 
relation lie on a line that passes outside of the origin. We show that we can remove all 
such holes from R to find an equivalent constraint language that is essentially convex, and 
thereby solve the problem in polynomial time. 

Combining these algorithmic results with matching NP-hardness results and the fact 
that CSP(r) is always in NP for a semilinear constraint language P (cf. Theorem 5.2 in 
Bodirsky et al. m) yields a dichotomy: 

Theorem 2. Let T be a finite set of semilinear eonstraints that eontains R^. Then, 
CSPfr) is either in P or NP-complete. 

Our result immediately generalises the first part of Theorem [TJ It also generalises 
another result by Bodirsky et al. [3] concerning expansions of {^2+} with relations that 
are first-order definable in {i2+}. One may note that this class of relations is a severely 
restricted subset of the semilinear relations since it admits quantifier elimination over 
the structure {-|-, {0}}, where -|- denotes the binary addition function. This follows from 
the more general fact that the first-order theory of torsion-free divisible abelian groups 
admits quantifier elimination (see e.g. Theorem 3.1.9 in [IB])- One may thus alternatively 
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view relations that are first-order definable in {-R+} as finite unions of sets defined by 
homogeneous linear systems of equations. 

We continue by generalising the second part of Theorem [U too: if T is semilinear and 
contains {i2+, {!}}, then the problem of optimising a linear polynomial over the solution 
set of CSP(r) is polynomial-time solvable if and only CSP(r) is polynomial-time solvable 
(and NP-hard otherwise). We also study the problem of finding integer solutions to CSP(r) 
for certain semilinear constraint languages P. Here, we obtain some partial results but a 
full classification remains elusive. Our results shed some light on the scalability property 
introduced by Jonsson and Loow m- 

This article has the following structure. We begin by formally defining constraint sat¬ 
isfaction problems and semilinear relations together with some terminology and minor 
results in Section [5J The algorithms and tractability results that are presented in Sec¬ 
tion O while the hardness results can be found in Section 01 By combining the results from 
Section [3] and m we prove Theorem [2] in Section [SI We partially generalise Theorem |2] to 
optimisation problems in Section jSl and we study the problem of finding integer solutions 
in Section [71 Finally, we discuss some obstacles to further generalisations in Section [51 
This article is a revised and extended version of a conference paper m- 

2. Preliminaries 

2.1. Constraint satisfaction problems. Let P = {i?i,..., Rn} be a finite set of finitary 
relations over some domain X (which will usually be infinite). We refer to T as a constraint 
language. In order to avoid some uninteresting trivial cases, we will assume that all 
constraint languages are non-empty and contain non-empty relations only. 

A first-order formula is called primitive positive if it is of the form 3xi,..., Xn-ipi A • • • A 
'0m.) where each 0j is an atomic formula, i.e., either x = y or i2(xq,... ,Xi^) with R the 
relation symbol for a A;-ary relation from P. We call such a formula a pp-formula. Note 
that all variables do not have to be existentially quantified; if they are, then we say that 
the formula is a sentence. Given a pp-formula $, we let Vars(‘h) denote the set of variables 
appearing in <h. The atomic formulas i?(xq ,... ,Xi^) in a pp-formula $ are also called the 
constraints of <h. 

The constraint satisfaction problem for a constraint language P (CSP(r) for short) is 
the following decision problem; 

Problem: CSP(P), where P is a finite set of relations over a domain X. 

Input: A primitive positive sentence over P. 

Output: ‘yes’ if is true in P, ‘no’ otherwise. 

The exact representation of the relations in P is unessential since we exclusively study 
finite constraint languages. 

A relation R{xi ,..., Xjf) is pp-definable from P if there exists a quantifier-free pp-formula 
ip over P such that 

R{xi, ...,Xk)^ 3yi,.. .,yn.p{xi,.. .,Xk,yi,.. .,yn)- 

The set of all relations that are pp-definable over P is denoted by (P). The following 
easy but important result explains the role of primitive positive definability for studying 
the computational complexity of CSPs. We will use it extensively in the sequel without 
making explicit references. 

Lemma 1 (Jeavons [IS])- Let F be a constraint language and P' C (P) a finite subset. 
Then CSP{F') is polynomial-time reducible to CSP{T). 

Let P = {i?i,..., Rk} and P' = {R'l ,..., R'k} be two constraint languages such that Ri 
and Rl are relations of the same arity. Given an instance <h of GSP(r), let denote the 
instance where each occurrence of a relation Ri is replaced by We say that GSP(r) 
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is equivalent to CSP(r') if is true in T if and only if <!>' is true in F'. It is clear that if 
CSP(r) and CSP(r') are equivalent CSPs, then they have the same complexity (up to a 
trivial linear-time transformation). 

2.2. Semilinear relations. The domain, X, of every relation in this article will be the 
set of rationals, Q, or the set of reals, M. In all cases, the set of coefficients, Y, will be 
the set of rationals, but in order to avoid confusion, we will still make this explicit in our 
notation. We define the following sets of relations. 

• LEx[Y] denotes the set of linear equalities over X with coefficients in Y. 

• LIx[Y] denotes the set of (strict and non-strict) linear inequalities over X with 
coefficients in Y . 

Sets defined by finite conjunctions of inequalities from LIx[Y] are called linear sets. 
The set of semilinear sets, SLx[Y], is defined to be the set of finite unions of linear sets. 
We will refer to S'L(q[Q] and S'L]k[Q] as semilinear relations over Q and M, respectively. 
One should be aware of the representation of objects in LEx\Y] and LIx\Y] compared 
to SLx[Y]. In LEx[Y] and LIx\Y], we view the equalities and inequalities as syntactic 
objects which we can use for building logical formulas. Now, recall the definition of a 
linear set: it is defined by a eonjunetion of inequalities. However, a linear set is not a 
logical formula, it is a subset of X^. The same thing holds for semilinear sets: they are 
defined by unions of linear sets and should thus not be viewed as logical formulas. This 
distinction has certain advantages when it comes to terminology and notation but it also 
emphasise a difference in the way we view and use these objects. The objects in LEx\Y] 
and LIx[Y] are often used in a logical context (such as pp-definitions) while the semilinear 
relations are typically used in a geometric context. 

Given a relation R of arity k, let R\x = R H X^ and r|x = {R\x | E T}. We 
demonstrate that CSP(r) and CSP(r|Q) are equivalent as constraint satisfaction problems 
whenever T C S'Lr[Q]. Thus, we will exclusively concentrate on relations from 5 'Lq[Q] in 
the sequel. Let T C 5L]r[Q] and let be an instance of CSP(r). Construct an instance 
of CSP(r|(Q) by replacing each occurrence of i? in ‘h by R\iq. If d>' has a solution, then 
has a solution since R\iq Y R for each i? E T. If <1> has a solution, then it has a rational 
solution by Lemma 3.7 in Bodirsky et al. [2] so has a solution, too. 

The following lemma is a direct consequence of our definitions: this particular property 
is often referred to as o-minimality in the literature In- 

Lemma 2. Let R E SLx\Y] be a unary semilinear relation. Then, R can be written as 
a finite union of open, half-open, and closed intervals with endpoints in Y U {—00,00} 
together with a finite set of points in Y. 

The set of semilinear relations can also be defined as those relations that are first- 
order definable in {!}} [E]- particular, SLx[Y] is closed under pp-definitions. 

Consequently, Lemma [2] is applicable to all relations discussed in this article. 

Lemma 3 (Lemma 4.3 in Bodirsky et al. [3]). Let ri,...,rfc,r E Q. The relation 
{(xi,...,Xfc) E I rixi -|- ... -|-TfcXfc = rj is pp-definable in {i?+,{l}} and it is pp- 
definable in {72+} ifr = 0. Furthermore, the pp-formulas that define the relations can he 
computed in polynomial time. 

It follows that L£'q[Q] C ({72+,{1}}) and L/q[Q] C ({72+, <,<, {!}}). One may also 
note that every homogeneous linear equation (with coefficients from Q) is pp-definable in 
{72+}. 

2.3. Unary semilinear relations. For a rational c, and a unary relation 1/ C Q, let 
c-[/ = {c-x|xEt/}E ({72+, [/}). When c = —1, we will also write —U for (—1) • U. 
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Given a relation R C and two distinct points a,b € we define 

C-R,a,b{y) = 3x1, R{xi, ...,Xk) A Af=i Xi = {l-y)-ai + y- bi. 

The relation Cji^a,b is a parameterisation of the intersection between the relation R and a 
line through the points a and b. Note that Cji^a,b is a member of (L£^q[Q] U {R}) so, by 
LemmaEl CR^a,b is a member of ({i?+, {!},-R}), too. 

A k-avy relation R is bounded if there exists an a € Q such that R C [—a, a]^. A unary 
relation U is unbounded in one direction if U is not bounded, but there exists an a € Q 
such that one of the following holds: U C [a,oo); or U C (— 00 , a]. A unary relation is 
called a bnu (for bounded, non-constant, and unary) if it is bounded and contains more 
than one point. 

Lemma 4. Let U be a unary relation in 5'Lq[Q] that is unbounded in one direction. Then, 

(1) ({R+, {!}, R}) contains a bnu. 

(2) if, in addition, U contains both positive and negative elements, then ({R+,R}) 
contains a non-empty bounded unary relation. 

Proof. (1) By Lemma El there exists an a > 0 such that either 

(1) (— 00 , —a] n R = 0 and [a, 00 ) C R; or 
(ii) (— 00 , —a] C R and [a, 00 ) fl R = 0. 

Assume that (0) holds. (The remaining case follows by considering —R.) By choosing 
a rational 5 > 2a, it is not hard to see that the relation 

R'(x) = 3y . y = b — X AU{x) AU{y) 

is bounded and contains an interval. The result then follows from Lemma El 

(2) Assume that (0) holds and let c € R be a negative element. (The remaining case 
follows by considering —R.) Then, 

U"{x) = 3y.ay = cxAU{x)A U{y) 

is bounded and contains the element a. The result again follows from Lemma El □ 

For a unary semilinear relation T C Q, and a rational 5 > 0, let T + X(5) denote the set 
of unary semilinear relations R such that T C U and for all x G R, there exists a y € T 
with |x — y| < 6. 

Example 1. The set {—1,1} +X(^) contains all unary relations U such that {—1,1} C 
Rc(-|,-i)u(l,|). 

Lemma 5. Let U ^ 0 be a bounded unary semilinear relation such that U fl (— 00 ,0) = 0. 
Then, {{R+, R}) contains a relation Us G {1} +X(5), for every rational 5 > 0. 

Proof. Let R+ = supR and U~ = infR. By Lemma El there exist elements p~^,p~ G R 
with R+ — < 6U^ and p~ — U~ < 6U~. The relation Us := {p~)~^ ■ U fl (p+)“^ • R 

is pp-definable in {R+,R} and satisfies: 1 G Us, supR^ < 1 + 6U~^{p'^)~^ < 1 + <5, and 
mfUs>l-5U-{p-)-^>l-5. □ 

Lemma 6. Let U be a bounded unary semilinear relation such that U n (—e, e) = 0 for 
some e > 0 and U fl —R A 0- Then, ({R+, R}) contains a relation Us G {—1,1} -\-T{5), 
for every rational (5 > 0. 

Proof. Let T = U f] —U. The proof then follows using a similar construction as in the 
proof of Lemma El D 
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2.4. Essential convexity. Let i? be a fe-ary relation over Q. The relation R is convex if 
for all p,q (z R, R contains all points on the line segment between p and q. We say that 
R is essentially convex if for all p,q & R there are only finitely many points on the line 
segment between p and q that are not in R. 

We say that R excludes an interval if there are p,q £ R and real numbers 0 < (5i < (^2 < 1 
such that p+ (q — p)y 0 R whenever 5i < y < 82 - Note that we can assume that (5i, <52 are 
rational numbers, since we can choose any two distinct rational numbers 71 < 72 between 
5i and 82 instead of 81 and < 52 . 

If R is not essentially convex, and if p and q are such that there are infinitely many 
points on the line segment between p and q that are not in R, then we say that p and q 
witness that R is not essentially convex. Due to Lemma [21 we conclude that a semilinear 
relation is essentially convex if and only if it does not exclude an interval. We say that a 
constraint language is essentially convex if all its relations are essentially convex. 

Theorem 3 (Theorem 5.1 and 5.4 in Bodirsky et al. [2|). IfT is a finite set of essentially 
convex semilinear relations, then CSPiT) is in P. 

3. Tractability 

In this section, we present our two main sources of tractability. Section 13.11 contains a 
new algorithm for semilinear constraint languages T containing {i?+,{l}} and such that 
(T) does not contain a bnu. In Section 13.21 we extend the applicability of Theorem |3] from 
essentially convex semilinear constraint languages to a certain class of semilinear CSPs 
that are not essentially convex. 

3.1. AfRne consistency. Instead of computing the exact solution set to a CSP instance, 
our approach will be to reduce an upper bound on this set as far as possible. In particular, 
we will maintain a representation of an affine subspace that is guaranteed to contain 
the solution set, and repeatedly intersect this subspace with every constraint in order 
to attempt to reduce it further. This can be seen as a form of local consistency. If we 
manage to reduce the upper bound to an empty set, then we are certain that the instance 
is unsatisfiable. We will show that under certain conditions, the converse holds; if the 
upper bound is non-empty, then there are necessarily solutions. To formalise this idea, we 
will need some definitions. 

For a subset S C Q”, let aff(S') denote the affine hull of S in Q”: 
k k 

aff(S') = (^XiPi \ k>l,Xi£ Q,pi £ S,'^Xi = 1}. 

i=l i=l 

An affine subspace is a subset S C Q" for which aff(5') = S. The points pi,... ,pk £ Q” 
are said to be affinely independent if xipi -!-••• + XkPk = 0 with xi + ■ ■ ■ + Xk = ^ implies 
xi = • • • = Xfc = 0. The dimension, dim(S'), of a set S' C Q"- is defined to be one less than 
the maximum number of affinely independent points in S. 

We define a notion of consistency for sets of semilinear constraints which we call affine 
consistency. Let P be a finite set of variables and let n = |P|. A set of constraints 
... ,Xif) with {x 7 ,... C P is affinely consistent with respect to a non-empty 

affine subspace 0 A C if aff(i?j r\A)=A for all i, where Ri := {(xi,..., Xn) £ | 

(x 7 ,... ,XiJ G Ri}. 

Algorithm [T] establishes affine consistency for a set of constraints and answers “yes” if 
the resulting affine subspace is non-empty and “no” otherwise. In the rest of this section, 
we show that this algorithm correctly solves CSP(r) when {i?+, {!}} C T is a semilinear 
constraint language such that (T) does not contain a bnu. Furthermore, we show that 
the algorithm can be implemented to run in polynomial time when applied to constraint 
languages of this kind. 
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Algorithm 1: Affine consistency 


Input: A set of constraints {Ri{xi^,... , 2 : 4 )} over variables V 

Output: “yes” if the resulting affine subspace is non-empty, “no” otherwise 


1 A := 

2 repeat 

3 foreach constraint Ri(xi^, ■ ■ ■, Xi^) do 

4 I A := aff(i?i n A) 

5 end 

6 until A does not change 

7 if A / 0 then return “yes” else return “no’ 


We begin by proving a technical lemma which is the basis for these results. 

Lemma 7. Let P = Pi U • • • U Pfc, Q = Qi U • • • U € 5Lq[Q] he two n-ary relations and 
Pi,, Pfc, Qi,... ,Qi linear sets. Assume that neither (LPq[Q]U{P}) nor (LPq[Q]U{(3}) 
contains a bnu. If aff(P) = aff((5) =: A, then aff(Pi (iQj) = A for some i and j. 

Proof. The proof is by induction on the dimension d = dim(A). For d = 0, both P and Q 
consist of a single point p. Clearly, Pj = {p} for some i and Qj = {p} for some j. Now 
assume that d > 0 and that the lemma holds for all P', Q' with aff(P') = sS.{Q') = A' 
and dim(A') < d. Let Po,Pi, ■ ■ ■ ,Pd be d + 1 affinely independent points in P and let 
go, gi, • • •, gd be d -|- 1 affinely independent points in Q. For 1 < i < d, consider the lines 
L? through po and pi, and the lines Lf through go and g^. Let H = {y € Q"' \ a ■ y = 0} 
{a G Q"') be a hyperplane in through the origin that is not parallel to any of the lines 
Lf or Lf. Then, H intersects each of the 2d lines. Let H{c) = {y G Q"' \ a - y = c} and let 
P(c) = {y G Q" I a • g 0 [-c, c]}. 

Express the line Lf as {y G Q” | y = {1 — x) ■ a + x ■ b, x G Q}, for some a,b G Q"". 
Let T = Cp^a,b € (LPq[Q] U {P}). Since T contains po and pi, it follows that T is not a 
constant and hence unbounded. By Lemma |U[^1), T is unbounded in both directions. By 
Lemma [21 P(cf) LI 'P T C P, for some positive constant cf. An analogous argument 
shows that that P(cJ) C Q, for some positive constant cj. Let c' be a positive constant 
such that po, go 0 B{c') and let c = max{c'} U {cf, Cj \ 1 < i,j < d}. This ensures that for 
any x > c, H{x)r\P intersects the lines Pf in d affinely independent points, and H (x) fl Q 
intersects the lines Pj in d affinely independent points. 

Let P' = H{x) n P, = H{x) n Pi, Q' = H{x) n Q, and Q' = H{x) n Qp We 
now have aff(P') = aff(Q') = A'(x) with dim(A'(x)) = dim(A) — 1, for every x > c. By 
induction on P' = P{ U • • • U P^ and Q' = Q[ U ■■■ U Q'l, it follows that aff(P(x) fl {Pi(x} C 
Qj{x})) = Qj{x)) ~ A'{x) for some i{x) and j{x). This holds for all x > c, 

hence there exist distinct xi,X 2 > c with i{xi) = i{x 2 ) = i' and j{xi) = j{x 2 ) = j'■ Since 
A'(xi), A'(x 2 ) C aff(Pi/ n Qj'), A'(xi) n A'(x 2 ) = 0, and dim(A'(x 2 )) = d — 1 > 0, it 
follows that aff(Pi/ flQp) strictly contains A'(xi), so we have A' C aff(Pi/ flQp) C A, and 
dim(A'(xi)) = dim(A) — 1. Therefore we have the equality aff(Pi/ n Qj/) = A. The lemma 
follows. □ 

For a semilinear relation R, we let size(P) denote the representation size of R, i.e., the 
number of bits needed to describe the arities and coefficients of each inequality in some 
fixed definition of R. 

Lemma 8. Let R G 5Pq[Q] be a relation such that (PPq[Q]U{P}) does not contain a bnu 
and let A C Q” be an affine subspace. Algorithmic computes a set of linear inequalities S 
defining aff(P fl A) in time polynomial in size(P) -|- size(A) and with size(S') < size(P) + 
size(A). 
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Algorithm 2: Calculate aff(i?n A) 

Input: A semilinear relation R = RiL) ■ ■ ■ U and an affine subspace A. 
Output: A set of inequalities defining aff(-R H A), or _L if aff(-R fl A) = 0. 

1 Find i that maximises d* := dim(aff(i2j n A)). 

2 if aS{Ri n A) = 0 then return _L 

3 Let I be the set of inequalities for Ri and J be the set of inequalities for A. 

4 S:=IUJ 

5 foreach inequality i € / U J do 

6 if dim(aff(S' \ {t})) = di then 

7 I 5:=5\{t} 

8 end 

9 end 

10 return S 


Proof. Let R = RiU ■ ■ ■ U R^ be the representation of R as the union of linear sets Ri. By 
Lemma [71 there exists an i such that aff(i? n A) = aff(-Rj H A) and since aS{Rj fl A) C 
aff(/Z n A) for all j, the algorithm will find such an i on line 1 by simply comparing the 
dimensions of these sets. If aff(i? n A) = 0, then the algorithm returns _L, signalling that 
the affine hull is empty. 

Otherwise, the affine hull of a non-empty polyhedron can always be obtained as a subset 
of its defining inequalities (cf. Schrijver [201 Section 8.2]). Here, some of the inequalities 
may be strict, but it is not hard to see that removing them does not change the affine hull. 
If 1 € / U J is an inequality that cannot be removed without increasing the dimension 
of the affine hull, then it is clear that i still cannot be removed after the loop. Hence, 
after the loop, no inequality in S can be removed without increasing the dimension of the 
affine hull. It follows that S itself defines an affine subspace, A 5 , and As = aff(A 5 ) = 
aff(i?inA) =aff(i?nA). 

Using the ellipsoid method, we can determine the dimension of the affine hull of a polyhe¬ 
dron defined by a system of linear inequalities in time polynomial in the representation size 
of the inequalities [201 Corollary 14.If]. To handle strict inequalities on line 1, we can per¬ 
turb these by a small amount, while keeping the representation sizes polynomial, to obtain 
a system of non-strict inequalities with the same affine hull. The algorithm does at most 
|7 U Jj -|- A: affine hull calculations. The total time is thus polynomial in size(7?) -|- size(A). 
Finally, the set 5 is a subset of / U J, so size(A 5 ') < size(i?) -|- size(A). □ 

Theorem 4. Let {i?+,{l}} C F C S'Lq[Q] he a finite constraint language. If there is no 
bnu in (F), then Algorithmic correctly solves CSPfT) and can be implemented to run in 
polynomial time. 

Proof. Assume that each relation 7? G F is given as 7? = 7?i U • • • U 7?^, where Ri is a linear 
set for each i. First, we show that the algorithm terminates with A equal to the affine 
hull of the solution space of the constraints. 

Assume that the input consists of the constraints 7?j(xi^,..., over variables V, 
i = 1,..., m. Let Z = denote the solution space of the instance. It is clear that Z 

is contained in A throughout the execution of the algorithm. Therefore, aff(Z) = aff(ZnA) 
so it suffices to show that aff(Z n A) = A. We will show that aff(ni=i .Ri H A) = A for 
all j = 1,... ,m. When the algorithm terminates, we have aff(Ri n A) = A for every 
i = 1,... ,m. In particular, the claim holds for j = 1. Now assume that the claim 
holds for j — 1. Then, P = Plti Rj C A and Q = Rj H A satisfy the requirements of 
Lemma |7| with aff(R) = aS{Q) = A. Therefore, we can use this lemma to conclude that 
aff (nLi Ri n A) = aff (R nQ) = A. 
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Finally, we show that the algorithm can be implemented to run in polynomial time. The 
call to Algorithm [2] in the inner loop is carried out at most mn times, where n = \V\. The 
size of R is at most size(i2) + logn, so the size of A never exceeds 0{mn{size{R) + logn)), 
where Ris a relation with maximal representation size. Therefore, each call to Algorithmic] 
takes polynomial time and consequently, the entire algorithm runs in polynomial time. □ 

3.2. Essential convexity. We will now identify another family of polynomial-time solv¬ 
able semilinear CSPs. This time, we base our result on essentially convex semilinear 
constraint languages (Theorem Ej) . We extend this result to the situation where we are 
only guaranteed that all unary relations that are pp-dehnable in the language are essen¬ 
tially convex. The idea is that even if we do not have the constant relation {1} to help 
us identify excluded intervals, we are still able to see excluded full-dimensional holes. We 
follow up this by showing that we can remove certain lower-dimensional holes and thus 
recover an equivalent essentially convex constraint language. We remind the reader that 
the dimension of a set is defined with respect to its affine hull, as in Section 13.11 

For € Q^, we let ||x|| denote the euclidean norm of x, and dist(x,y) = ||x — y|| the 
euclidean distance between x and y. 

Lemma 9. Let U € {1} + T(c) for some 0 < c < 1 and assume that R € S'Lq[Q] is a 
semilinear relation such that every unary relation in {{Rj^,U,R\) is essentially convex. 
Then, R can he defined by a formula ipo A A ■■■ A ^ipk, where ipo defines a convex 
semilinear set, and ipi,...,ipk ore conjunctions over L/q[Q] that define convex sets of 
dimensions strictly lower than the dimension of the set defined by ipo . 

Proof. Let conv(i2) denote the convex hull of R and let d denote its dimension. The 
set conv(i?) is semilinear (see, for instance, Stengle et al. [21ji. Let (po be a formula for 
conv(i?) and let yji V • • • V be a formula for conv(i?) \ i? on quantifier-free DNF over 
T/q[Q]. It remains to show that for each i, the dimension of the convex set Si defined by 
ifi is smaller than d. To prove this, we show that for every point p in Si, and every e > 0, 
there exists a point x in R such that dist(p, x) < e. Since every d-dimensional convex set 
contains a small d-dimensional open ball around every point in its interior, it follows from 
this that none of the sets Si can be d-dimensional. 


Xm 



Figure 1. An illustration of the entities involved in the induction step. 

Caratheodory’s theorem (cf. Schrijver [501 Section 7.7]) states that for every p G conv(i?), 
we can find m -|- 1 < d -|- 1 affinely independent points, xq, ..., Xm € R, such that p lies in 
B = conv({xo,... ,Xm})- By induction over m, we show that for every point b £ B, and 
every e > 0, there is a point z £ R such that dist(6, z) < e. For m = 0, this statement 
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follows trivially as each Xj was chosen from R. Now assume that 0 < m < d and that the 
statement holds for all 0 < m' < m. By the induction hypothesis, the statement holds for 
the set A = conv({xo, ■ ■ ■, Xm-i})- Every b £ B can be written as b = yi, ■ a + {1 — yi,) ■ Xm 
for some a £ A and 0 < 7/fe < 1. Let a' £ Rhe a point in A that is at distance at most e/2 
from a and let 6' = y;, • a' + (1 — y;,) • Xm- Then, 

dist(6,6') = II( 2/6 ■ o + (1 - yh) ■ Xm) - {Vh • o' + (1 - y^) ■ Xm)\\ 

< yfello - a'll 

< e/2. 

Let 5 > 0 be a small constant to be fixed later. By Lemma El there exists a unary 
relation Us £ {1} +Z{5) fl (i?+, U). Consider the following relation: 

r(y) = 3t, z.Usit) A R{z) A z = y ■ a' + {t - y) ■ Xm- 

Since Us £ ({i2+,C/}), we also have T £ ({i?+,[/}). By assumption, T does not exclude 
an interval, so there exists a y^ such that T{y'^) and |y^ — y;,! < 6 . Then, by the definition 
of T, there exists a t G (1 — 5,1 + (5) and a point z £ R such that: 

dist(6', z) = II (yb • o' + (1 - ys) ■ Xm) - {y'b ■ a' + {t - y'^) ■ Xm) || 

= WiVb - y'b) ■ a' + {1-t) ■ Xm + {y's - yb) ■ Xm\\ 

< 11(2/6 - y'b) ■ a'll + 11(1 - i) ■ Xm\\ + \\{y'b - yb) ■ Xm\\ 

< {\yh - y'b\ + |1 - + \y'b - yb\) niax{||a'||, ||xm||} 

< 36C, 

where C := max{||a'||, ||xm||, 1} is a constant for a fixed B, and the first inequality follows 
from the triangle inequality. The claim now follows for the point b by taking <5 = (e/2) • 
(3(7)“^ since dist(5, < dist( 6 ,5') + dist( 6 ', 2 ;) < e. □ 

Theorem 5. Let {ii+} C L C 5 Lq[Q] be a constraint language. Assume that there exists 
a unary relation U £ {1} +T(c) fl (T), for some 0 < c < 1, and that every unary relation 
in (r) is essentially convex. Then, CSPiT) is equivalent to CSP{T') for an essentially 
convex constraint language T' C S'LqIQ]. 

Proof. If T is essentially convex, then there is nothing to prove. Assume therefore that 
r is not essentially convex. By Lemma O each R £ T can be dehned by a formula 
(fo A A ■ ■ ■ A ^ipk, where ipQ, (fi, ..., (pk are conjunction over L/q[Q], and (fi, ..., pk 
define sets whose affine hulls are of dimensions strictly lower than that of the set defined 
by pq. Assume additionally that the formulas are numbered so that the affine hulls of the 
sets dehned by pi,... Pm do not contain ( 0 ,... , 0 ) and that the affine hulls of the sets 
dehned by Pm+i, ■ ■ ■ ,Pk do contain (0,...,0). Dehne R' by 

p A ^p'l A • • • A ^p'.^ A ^Pm+l A • • • A ^pk, 

where p[ dehnes the affine hull of the set dehned by pi. Then, the constraint language 
r' = {R' I i? G r} is essentially convex since witnesses of an excluded interval only occur 
inside an affine subspace not containing ( 0 ,... , 0 ); otherwise we could use such a witness 
to pp-dehne a unary relation excluding an interval. 

Let be an arbitrary instance of CSP(r) over the variables V = {xi,... ,Xn} and 
assume $ = 3xi,..., where "f/ is quantiher-free. Construct an instance <!>' of CSP(r') 
by replacing each occurrence of a relation R in by R'. Clearly, if is satishable, then so 
is <h. Conversely, let s £ be a solution to ‘b and assume that ‘b' is not satishable. Let L 
be the line in through ( 0 ,..., 0 ) and s and let U be the unary relation ^ (T). 

All tuples in U correspond to solutions of that are not solutions to <b'. 

Fix a constraint R{xi,... ,xi) in $ and consider the points in U that satisfy this con¬ 
straint but not R'{xi,... ,xi). These are the points p £ on L for which (p(xi),... ,p{xi)) 
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satisfies ((^i V- • •V(/3^) A-'(<y9iV- • -Mipm)- For each 1 < i < m, ip'^ satisfies at most one point 
on L since otherwise the affine hull of the relation defined by pi would contain (0,..., 0). 
Hence, each constraint in d> can account for at most a finite number of points in U, so U 
is finite. 

Assume hrst that \U\ > 1. Then, U is not essentially convex which contradicts the 
assumption that every unary relation in (F) is essentially convex. Assume instead that 
U = {!}, where the single point in U corresponds to the solution s. Recall that F is not 
essentially convex. Let R E F be a k-aiy relation that is not essentially convex and let 
p,g E witness this. The relation E (F) since {1} E (F). Then, E (F) is 

unary and not essentially convex which leads to a contradiction. It follows that if ‘h is 
satishable, then so is ‘hh □ 


4. NP-hardness 

We now derive a unified condition for all hard CSPs classified in this article. It is based 
on a polynomial-time reduction from the NP-hard problem Not-All-Equal 3SAT m, i-e. 
the problem CSP({Rnae}) where Rnae = {~1) 1}^ \ {(~1) ~1) ~1)) (Ij 1; !)}• The proof 
is divided into three different lemmas. First, we present a reduction from Not-All-Equal 
3SAT to a simple semilinear CSP. We then show that having a bnu T that is bounded 
away from 0 allows us to pp-define unary relations that are, in a certain sense, close to 
being either the relation {1} or {—1,1}. In the final step, we combine these two results 
and show that having a bnu T that excludes an interval and that is bounded away from 
0 is a sufficient condition for CSP({R+,T}) to be NP-hard. 

Lemma 10. Let T E {—1,1} -|-X(^). Then, CSP{{Rj,.,T}) is NP-hard. 

Proof. The proof is by a polynomial time reduction from CSP({Rnae})- Let 4> denote 
an arbitrary instance of CSP({Rnae})- Construct an instance <h' of CSP({R+,T}) as 
follows. Impose the constraint T{v) on each variable. Eor each constraint z) in 

<h, introduce the constraints x-\-y + z-\-w = 0 and T{w), where tc is a fresh variable. 

Assume that has a solution. Consider a constraint R^p,:^{x,y, z) in <I>. If two of the 
variables are assigned the value 1, then the equation x-|-y-|- 2 ;-|-t(; = 0is satished by 
choosing w = —1. If two of the variables are assigned the value —1, then the equation 
x-|-y-|- 2 :-|-tc = 0is satished by choosing w = 1. Hence, <h' is satishable. 

Assume that has a solution s'. Then, <I> has a solution s dehned by s{x) = 1 if 
s'(x) > 0 and s{x) = —1 if s'{x) < 0. Assume to the contrary that s(x) = s(y) = s(z) = 1 
for some variables with a constraint R]s;ae(x, y, z). Consider the equation x + y-\-z + w = 0 
in <h'. By the assumption on T, we have s'{x) s'{y) s'{z) > |, and hence s'{w) < —|. 

But this is a contradiction as the constraint T{w) is also in <!>'. We can similarly rule out 
the case s(a:) = s{y) = s{z) = —1. This proves that s is a solution to 4>. □ 

Lemma 11. Let T ^ 0 be a bounded unary relation such that T n (—e,e) = 0 , for some 
e > 0. Then, either (R_|_,T) contains a unary relation Us E {1} for every 5 > 0; or 

(R+,T) contains a unary relation Us E {—1,1} +T{5), for every <5 > 0. 

Proof. If r n —T 7 ^ 0 , then the result follows from Lemma EJ Otherwise, by Lemma [21 
there exists a constant c"*“ > 0 such that the set T"*" = {x E T | |x| > c"*“} is non-empty 
and contains points that are either all positive or all negative. Similarly, there exists a 
constant c“ > 0 such that T~ = {x E T | |x| < c~} is non-empty and contains points 
that are either all positive or all negative. Let a E T'*' and b E T~. Assume that both 
sets contain positive points only or that both sets contain negative points only. Then, the 
result follows using LemmaElwith the relation U = a~^ ■ T n b~^ ■ T (or —U if the points of 
U are negative). The case when the one set contains positive points and the other contains 
negative points is handled similarly using the relation U' = a~^ ■ T n h~^ ■ {—T). □ 
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Lemma 12. Let T be a bnu such that T n (—e, e) = 0 , for some e > 0, and U be a unary 
relation that excludes an interval. Then, CSP{{Rj^,T,U}) is NP-hard. 

Proof. We show that T, U) contains a unary relation {—1,1} +X(^). The result then 
follows from Lemma fTUl If already (i?+,T) contains such a relation, then we are done. 
Otherwise, by Lemma fTTl (ii+,T) contains a unary relation Us G {1} +T{5), for every 
h > 0. Since U excludes an interval, there are points p,q ^ U and 0 < hi < ^2 < 1 such 
that p+ (q — p)y ^ U whenever hi < y < h 2 . Furthermore, p and q can be chosen so that 
hi < 1/2 < h 2 , and by scaling U, we may assume that \q — p\ = 2. Let m = {p + q)/‘i. 
Note that T n (m — e',m + e') = 0, for some e' > 0. Similarly, possibly by first scaling T, 
let p', q' € T he distinct points with \q' — p'\ = 2 and let m' = {p' + q')/2. 

Now, define the following unary relations: 

To{x) = 3y3z . Us{y) A z = x — y ■ m A U(z) 

Too{x) = 3y'3z'. Us{y') A z' = x - y' ■ m' A T{z'). 

The relations Tq and Tqo are roughly translations of U and T, where the constant relation 
{1} has been approximated by the relation Us. Since 1 € Us, we have {—1,1} C To,Too. 
Hence, if h is chosen small enough, then the relation Tq H Tqo € (i?+, T, U) will satisfy the 
conditions of Lemma El This finishes the proof. □ 

5. Semilinear expansions of { R +} 

In this section, we prove our main result: Theorem El We divide the proof into two 
parts. Consider the following two properties: 

(Pq) There is a unary relation U in (T) that contains a positive point and satisfies 
[/ n (0, e) = 0 for some e > 0. 

(Poo) There is a unary relation U in (P) that contains a positive point and satisfies 
U n (M, oo) = 0 for some M < oo. 

In the first part of the proof (Section 15.1|) . we consider constraint languages that si¬ 
multaneously satisfy the properties (Pq) and (Poo)- In the second part (Section 15.2|) . we 
consider constraint languages that violates at least one of them. In both parts, we give 
a detailed description of the boundary between easy and hard problems. By combining 
Theorem El and Theorem El we establish Theorem El 

In addition to the two algorithmic results in Section KLll there is also a trivial source 
of tractability. A relation is 0-valid if it contains the tuple (0,..., 0) and a constraint 
language is 0-valid if every relation in it is 0-valid. Every instance of a CSP over a 0-valid 
constraint language admits the solution that assigns 0 to every variable. 

When we consider constraint languages that are not 0-valid, the following lemma shows 
that there is always a pp-definable unary relation that is not 0-valid. 

Lemma 13. Let {12+} CPC S'Lq[Q] be a constraint language. IfT is not 0-valid, then 
(P) contains a non-empty unary relation that is not 0-valid. 

Proof. By assumption, P contains some k-aiy relation R that is not 0-valid, and by our 
definition of a constraint language, R is non-empty. Let t G 12 be a tuple that contains 
the largest possible number m of zeroes. Assume for simplicity that the first m entries of 
t equals 0. Consider the following unary relation in (P). 

1/ = (x G Q I 3ym+i ■ --yk-i ■ 12(0,0,.. .,0,ym+i, ■ ■ ■,yk-i,x)} 

The relation U is non-empty and not 0-valid. □ 
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5.1. The case (Pq) and (Pqo)- The following theorem covers the case when the con¬ 
straint language satisfies both of the properties (Pq) and (Poo)- As a corollary, we obtain 
a complete classification for semilinear constraint languages containing {i?+,{l}}. The 
latter result is interesting in itself and it will also be used in Section 15.21 and Section El 

Theorem 6. Let {12+} CPC S'Lq[Q] be a finite constraint language that satisfies (Pq) 
and (Poo)- The problem CSP{T) is in P if 

• r is 0-valid (trivially); 

• (P) does not contain a bnu (by establishing affine consistency); or 

• all unary relations in (P) are essentially convex (by a reduction to an essentially 
convex constraint language). 

Otherwise, CSP{T) is NP-hard. 

Proof. Let U be the set of all bounded, non-empty unary relations U in (P) such that 
U n (—e,e) = 0 for some e > 0. Assume that P is not 0-valid. First, we show that U is 
non-empty. By Lemma fTHl (P) contains a non-empty unary relation that is not 0-valid. 
Scale this relation so that it contains 1 and call the resulting relation U'. Let Uq G (P) 
be a unary relation witnessing (Pq) and let Uoo G (P) be a unary relation witnessing 
(Poo)- Scale Uq and Uoo so that some positive point from each coincides with 1 and let 
T = U' n Uq n Uoo- If T does not contain a negative point, then T £ U. Otherwise, T 
contains a negative point b. It follows that T n 6 - T £U. Hence, the set IL is non-empty. 

Assume that (P) does not contain a bnu. Then, neither does U and hence U contains 
only constants. It follows by Theorem 0] that establishing affine consistency solves CSP(r). 

Otherwise, U contains a bnu. If all unary relations of (P) are essentially convex, then 
by Lemma [TT] and Theorem El CSP(r) is equivalent to CSP(r') for an essentially convex 
constraint language Tractability follows from Theorem [3l 

Finally, if U contains a bnu and (P) contains a unary relation that excludes an interval, 
then NP-hardness follows from Lemma fT^ □ 

Corollary 1. Let {i2+,{l}} CPC S'LqIQ] be a finite constraint language. The problem 
CSP(T) is in P if (P) does not contain a bnu or if P is essentially convex. Otherwise, 
CSP{T) is NP-hard. 

Proof. If (P) does not contain a bnu, then tractability follows from Theorem [H If all 
relations in P are essentially convex, then tractability follows from Theorem EJ 

Otherwise, (P) contains a bnu, and since {i2+,{l}} C P, (P) also contains a unary 
relation that is not essentially convex. Since {1} G P is not 0-valid, NP-hardness then 
follows from Theorem El □ 

5.2. The case -'(Pq) or -i(P(X)). Let {22+} CPC S'Lq[Q] be a constraint language such 

that either (Pq) or (Poo) is violated. In this section, we show that P can be replaced by 
an equivalent constraint language of a restricted type. Let 225 Lq[Q] denote the set of 
relations that are are finite unions of homogeneous linear sets. We will call such relations 
homogeneous semilinear relations. We remind the reader that we can always pp-define the 
relations {0} and M = {{x,—x) \ x G Q} in P: x = 0 22+(x,x,x) and {x,y) £ M 4^ 

i2+(x,y,0) 3z.R+{x,y, z) A R+{z, z, z). Hence, we can freely use the constant 0 and 

negation in forthcoming pp-definitions. 

From now on, let Q+ = {a G Q | a > 0}, Q_ = {a G Q | o < 0}, and Q^o = Q- U Q+ = 
Q\{0}. For a relation R £ 5Lq[Q], define cone(i2) = {A-x | A G Q+, x G 22} to be the cone 
over 22. For a constraint language P C S'Lq[Q], let cone(P) = {cone(22) | 22 G P}. Note 
that, for {22+} C P C 5Lq[Q], we have cone(P) C 22S'Lq[Q], and since cone(22+) = 22+, 
we also have 22+ G cone(P). 

For an assignment s : P —>■ Q and a rational c G Q, let c - s denote the assignment 
X !->■ c - s(x). 
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Theorem 7. Let {i?+} C r C S'Lq[Q] be a constraint language such that either (Pq) or 
(Poo) is violated. Then, CSP{T) is equivalent to CSP{cone{T)). 

Proof. Assume that T does not satisfy (Pq)- The proof for the case when P does not satisfy 
(Poo) follows similarly. 

Let i? be a relation in P and let <y9 = <y9i V • • • V be a quantifier-free DNF formula 
for R, where each formula (pj is conjunction of strict and non-strict inequalities. Remove 
every disjunct (pj that contains a non-homogeneous inequality which is not satisfied by the 
(0,..., 0)-tuple. Let S be the relation defined by the resulting formula (p' = p'lV ■ ■ ■ V p'^,. 
Since P does not satisfy (Pq), it follows that for every point x in R\S, there is a point x' in S 
that lies on the open line segment between (0,..., 0) and x. Therefore, cone(S') = cone(i?). 
Next, for each j, let Sj be the relation defined by p'j. Remove every non-homogeneous 
inequality from pb, let p'- be the resulting formula and let Tj be the relation defined by 
p'-. Clearly, cone(S'j) C cone(Tj). Let A-x be a point in cone(rj) with A G Q+ and x G Tj. 
Since every non-homogeneous inequality in p'j is satisfied by the (0,..., 0)-tuple, it follows 
that they are satisfied by every point in a small ball B centred at (0,..., 0). Let x' be a 
point in B on the line segment between (0,... , 0) and x and note that every homogeneous 
inequality in p'j satisfies x and therefore also x'. It follows that x' is in Sj so x and A • x 
are in cone(Sj), which shows that cone(rj) C cone(S'j). Let p" = y?" V • • • V p'y and let T 
be the relation defined by p". Then, cone(R) = cone(T) and cone(r) = T since p" only 
contains homogeneous inequalities. Therefore, p" defines cone(i?), so cone(R) G 
and P' C R'SLq[(Q]. 

Por the equivalence of CSP(P) and CSP(r'), arbitrarily choose an instance d> of CSP(r). 
Construct an instance <!>' of CSP(r') by replacing each occurrence of a relation R in <I> by 
cone(R). Every solution to is also a solution to <!>'. It remains to show that if <!>' has a 
solution, then so does 

Let s : Vars(d>') —>■ Q be a solution to d>' and assume without loss of generality that s 
is integral. If s = 0, then it follows immediately that s is a solution to since, for every 
R G P, (0,... , 0) G cone(R) if and only if (0,..., 0) G R. Assume therefore that s ^ 0. 
Let Ri(xi),... ,Rm(xm) be the constraints of <I>. For every j, s{xj) G cone(Rj) holds. By 
the construction of cone(Rj), this implies that r • s{xj) G Rj, for some r > 0. Define the 
unary relation U G (P) by the pp-formula V’(y) = . x = y • s{xj) A Rj{x). Now r G U, 

so by the assumption on P, it follows that (0, e^) C U, for some Sj > 0, and hence that 
y ■ s{xj) G Rj, for all y G (0, £j). Let e = min^ Sj. Then (e/2) • s is a solution to <I>. □ 

By Theorem[71 it is thus sufficient to determine the computational complexity of CSP(P) 
for {R+} C P C HSLq[Q]. 

Given a relation R C Q^, we say that a function e : Q —>■ Q is an endomorphism of R 
if for every tuple (oi,..., Ok) G R, the tuple (e(ai),..., e(afc)) G R. One may equivalently 
view an endomorphism as a homomorphism from R to R. We extend this notion to 
constraint languages P = {Ri,..., Rn}: a function e : Q —)• Q is an endomorphism of P if 
e is an endomorphism of Rj, 1 < i < re. 

Lemma 14. Let a > 0 be a rational number. Every R G has the endomorphism 

e(x) = a ■ X. 

Proof. We know that R can be written as R = Uj^i where Rj, 1 < i < rre, is defined 
by a (finite) system of homogeneous linear (strict or non-strict) inequalities. Consider an 
inequality Cj • Xj > 0 in such a system. We immediately see that 

n n n n 

y/ Cj • Xj > 0 a • Cj • Xj > 0 a ■ Ci ■ Xi > O Cj • e(xj) > 0. 

7=1 7=1 7=1 7=1 

This equivalence also holds if we consider strict inequalities. Therefore, each Rj, 1 < i < rre, 
has the endomorphism e. 
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Now, arbitrarily choose a tuple t = {ti,... ,tk) € R and assume that t € Hi. It follows 
that (e(ti),..., e(tfc)) € Hi C R, so the function e is an endomorphism of R. □ 

A direct consequence of Lemma fT^ is the following: if an instance of CSP(Lf5'LQ[Q]) 
has a solution s, then a • s is a solution for every rational number a > 0 . 

The complexity classification of constraint languages that violate either (Pq) or (Poo), 
in Theorem [SJ follows from two intermediate results which we now present in Lemma fTCl 
and Lemma m 

Lemma 15. Let T be a subset of HSLq[Q] and let U be a unary relation in (T). If U 
contains an element p > 0, then Q+ CP. If U contains an element p < 0, then Q_ C U. 

Proof. Let (7 G Q be any element with the same sign as p. By Lemma [TU e(x) = {q/p) ■ x 
is an endomorphism of U. Since p € U, it follows that q = e{p) £U. □ 

Lemma 16. Let {P+} C T C PS'Lq[Q] be a finite constraint language. Then, either 

• r is 0-valid; or 

• C'S'P(r) is polynomial-time equivalent to PS'P(r U {{1}}). 

Proof. Assume that T is not 0-valid. By Lemma fTHl (T) contains a non-empty unary 
relation that is not 0-valid. We consider three different cases. 

Case 1. (r) contains a non-empty unary relation P such that 0 0 P and P C Q+. By 
Lemma fl^ Q+ C P so P = Q+. We claim that CSP(r U {{1},Q+}) is polynomial-time 
equivalent to CSP(r U {Q+}). The polynomial-time reduction from right to left is trivial. 
To show the other direction, let d> be an arbitrary instance of CSP(ru{{l}, Q+}). Assume 
without loss of generality that the relation { 1 } appears in exactly one constraint {l}(x). 
Construct by replacing this constraint with Q_|_(x). 

If has no solution, then <1> has no solution. Suppose instead that has the solution 
s. Then we know that s(x) > 0. Choose o G Q such that a ■ s(x) = 1. By Lemma flTl the 
function o • s is then a solution to <I>. 

Case 2. (T) contains a non-empty unary relation P such that 0 0 P and P C Q_. By 
Lemma [El Q- C P so P = Q_. We can now pp-define Q+ since x > 0 —x < 0 and go 
back to Case 1. 

Case 3. (P) contains a non-empty unary relation P such that 0 ^ P and no unary relation 
P' G (P) equals Q+ or Q_. Lemma [El implies that P = Q_ U Q+. 

We claim that CSP(r) is polynomial-time equivalent to CSP(ru {{!}}). The reduction 
from left to right is trivial. To show the other direction, let = 3xi,... ,Xm-p{xi,... ,Xm) 
be an arbitrary instance of CSP(r U {{1},P}), where p is quantifier-free, and assume 
without loss of generality that the relation { 1 } appears in exactly one constraint {l}(xm). 
Construct by replacing this constraint with Q^o(3^m)- 

If has no solution, then ‘k has no solution. Suppose instead that has a solution. 
Assume first that every solution assigns a negative number to the variable Xm- Then we 
can pp-define a unary relation T C Q_ by 

T{Xm) = 3 X 1 , ... , Xm-l ■ p{xi, . . . , Xm) 

and this contradicts our initial assumptions. Thus, there is a solution s such that s{xm) > 
0. Choose a G Q such that a ■ s{x) = 1. By Lemma fTTl the function a ■ s is a solution to 

d>. □ 

Theorem 8 . Let {^2+} PTC 5Lq[Q] be a finite constraint language that violates (Pq) 
and/or (Poo)- The problem CSP{I') is in P if 

• r is 0-valid; 

• (cone(r) U {{!}}) does not contain a bnu; or 

• cone(r) is essentially convex. 
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Otherwise, CSPiT) is NP-hard. 

Proof. By Theorem[71 CSP(r) is equivalent to CSP(cone(r)). By LemmafTHl CSP(cone(r)) 
is either trivially in P, if it is 0-valid, or CSP(cone(r)) is polynomial-time equivalent to 
CSP(cone(r) U {{!}})• In the latter case, the result follows from Corollary[TJ □ 

6. Optimisation 

In this section, we study the optimisation problem where the objective is to maximise 
a linear function over the solution set of a semilinear CSP. For an arbitrary constraint 
language P C 5Lq[Q], we formally define the problem Opt(r) as follows. 


Problem: Opt(r) 

Input: A CSP(r) -instance and a vector c G QVars('i') _ 

Output: One of the following four answers. 

• ‘unbounded’ if for every K £ Q, there exists a solution x such that c^x > 

K. 

• ‘optimum: K' if there exists a AT G Q and a solution x such that c^x = K, 
but there is no solution x' such that c^x' > K. 

• ‘optimum is arbitrarily close to AT’ if there exists a AT G Q such that there 
is no solution x satisfying c^x > K, but for every K' < K there is a 
solution x' with c^x' > K'. 

• ‘unsatisfiable’ if there is no solution. 


By Lemma El the problem Opt({A+, <, {!}) is polynomial-time equivalent to linear 
programming. Bodirsky et al. [2] have shown that for semilinear constraint languages 
containing {A+, <,{!}}, the problem CSP(r) is polynomial-time solvable (NP-hard) if 
and only if the problem Opt(r) is polynomial-time solvable (NP-hard) (cf. Theorem [T]). 

In Theorem IIUI we show that, for semilinear constraint languages containing {A+, {!}}, 
the complexity of the decision problem and of the optimisation problem is similarly related. 
We first prove an analogue of Theorem 0] for the optimisation problem. 

Theorem 9. Let {A+, {!}} CPC SAQjQ] he a finite constraint language. If there is no 
bnu in (P), then Opt{T) can be solved in polynomial time. 

Proof. Let be an instance of CSP(r), let V = Vars(<I>) = {xi,... ,Xm}, and let c G 
be a vector. Assume = 3xi,... ,Xm ■ ‘P where ip is quantifier-free. Algorithm [T] in 
Section El finds the affine hull A of the set of satisfying assignments to ‘h in polynomial 
time. If A = 0 , then we answer ‘unsatisfiable’. 

Otherwise, the affine hull A is represented by a set of inequalities, each with repre¬ 
sentation size that is polynomial in the input size. Therefore, we can solve the system 
zi,Z 2 G A, c^{zi — Z 2 ) > 0, in polynomial time. Assume that this system has a solution. 
Let k = dim(A) -|- 1 and let yi,... ,yk be affinely independent satisfying assignments to 
$. Then, we can write zi = Y.i=i o-uUi and Z 2 = Y.i=i with “i* = J2i=i “ 2 * = 1- 
Since 

k k 

c^{zi - Z 2 ) = aiiC^yi - a 2 iC^yi > 0, 
i=l i=l 

we must have c'^yi 7 ^ c'^yj for some 1 < i,j < k. Let U = G (P), where 

= {(xi, ... ,Xn) G Q'^ I ifixi,... ,Xn) is true in P} and for each a € U, let ya G 
denote the corresponding point on the line through y^ and yj. Fix an arbitrary constant 
AT G Q. Since there is no bnu in (P), it follows from Lemma 01(1) that there is a point 
a £ U such that ya £ satisfies c'^ya > K. Since ya is a satisfying assignment, we can 
therefore answer ‘unbounded’. 
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Otherwise, c^{zi — Z2) = 0 for all zi, Z2 € A, so is constant for z ^ A. In polynomial 
time, we can find a z G A with polynomial representation size. It then suffices to evaluate 
c'^z and answer ‘optimum: c'^z’. □ 

Theorem 10. Let C F C S'Lq[Q] be a finite constraint language. The problem 

Opt(T) is polynomial-time solvable if (F) does not contain a bnu or if F is essentially 
convex. Otherwise, OptiT) is NP-hard. 

Proof. The polynomial-time solvable cases follow from Theorem [T] and Theorem O The 
hardness follows from Corollary [TJ □ 

A comparison between Theorem [TU] and Corollary [T] shows that, for a semilinear con¬ 
straint language F containing {i2+,{l}}, CSP(F) is polynomial-time solvable (NP-hard) if 
and only if Opt(F) is polynomial-time solvable (NP-hard). The following example shows 
that this tight relationship between the complexity of a constraint satisfaction problem 
and its corresponding optimisation problem cannot be further extended to the class of all 
semilinear constraint languages containing the relation ii+. 

Example 2. Let R = {(0,0,0,0)} U {{x,y,z,l) \ {x,y,z) G i?NAE} (cf. Section\^. Note 
that F = (i?, i?+} is semilinear, 0-valid, and that F satisfies both (Pq) and (Poo)- Let <I> be 
an arbitrary instance of CSP{{R^ ae})- Construct an instance of Opt{T) by introducing 
an auxiliary variable w, and for each constraint R^AEix,y, z) in <I>, introduce a constraint 
R{x,y,z,w) in d>'. Finally, let the vector c € QVars('J>) defined by c^ = 1 and Cx = 0 
for all other variables x. Then, the instance has a solution if and only if an optimal 
solution of has value 1. We conclude that CSP{T) is polynomial-time solvable (since F 
is 0-valid), but that Opt{T) is NP-hard. 

7. Integer solutions 

In this section, we study the problem of finding integer solutions to CSPs defined 
by semilinear relations. We consider two different approaches: (1) allowing an additional 
unary constraint that forces a chosen variable to take an integral value, and (2) identifying 
constraint languages which guarantee the existence of integer solutions. 

The reader should note that in the first approach we do not consider semilinear relations 
defined over the integers. Instead, we consider ways of checking whether a given problem 
instance has a solution where some variables are assigned integral values. Some of the 
problems in the second approach can be seen as semilinear CSPs over the integers, but 
our methods do not lend themselves to a systematic study of such CSPs. See [S] for a 
recent approach to such a systematic study. 

7.1. The relation Z. The unary relation Z can be used to ensure that a variable is given 
an integral value. By Lemma O this relation is not semilinear over Q, so the constraint 
languages that we classify in the next theorem are formally not semilinear. 

Theorem 11. Let {i?+} C F C be a finite constraint language that satisfies (Pq) 

and (Poo)- The problem CSP{T U {Z}) is in P if 

• F is 0-valid; or 

• (F) does not contain a bnu. 

Otherwise, CSP{T U {Z}) is NP-hard. 

Proof. If F is 0-valid, then F U {Z} is 0-valid, so every instance admits the solution 

( 0 , 0 ,..., 0 ). 

Otherwise, assume that (F) does not contain a bnu. Let be an arbitrary instance of 
CSP(F U {Z}), let I C Vars(<h) be the set of variables that are constrained by Z in <1>, and 
let d>' be the instance of CSP(F) obtained from <I> by removing all Z-constraints. By an 
argument on the set of all bounded, non-empty unary relations in (F) similar to that used 
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in the proof of Theorem [ 6 l it follows that (T) contains the relation {!}. Let S be the set 
of satisfying assignments to By running Algorithm [U we obtain a system of inequalities 
that defines the affine hull A of the satisfying assignments S. 

We now substitute each such inequality for an equality. The resulting system of linear 
equalities still defines A. Let A' = | x € A}. We can compute a system of linear 

equations for A' in polynomial time by first computing a parameter form for A, removing 
the coordinates not corresponding to I, and then computing the equivalent system of linear 
equations. This can be in polynomial time by being careful with the representation sizes of 
the intermediary results (cf. Schrijver [2D1 Section 3]). We then solve the resulting system of 
linear equations for an integer solution in polynomial time (cf. Schrijver [ 2 Ul Corollary 5.3]). 
If no such solution exists, then is unsatisfiable. Otherwise, the integer points in A' are 
given by L = {cq + I ■^ij • • • S Z}, for some linearly independent vectors 

Co,..., Cfc E Z^, where k = dim(A'). The vectors c* can be found explicitly in polynomial 
time, but since we are only interested in showing that there exists a satisfying assignment 
to ‘h, it suffices that L has the aforementioned form. 

For p ^ A' and constant e > 0, define B{p,e) = {x £ A' \ \\p — x|| < e}. Let 
S' = {niix) I X E S'} be the projection of S to the variables in I and note that S' E (T). 
By Lemma [71 S' contains a linear set R C such that aff(i?) = A'. Let p £ R and e > 0 
be such that B{p,£) C i? C S'. We claim that there exist distinct qi,q 2 £ L such that 
the line through qi and q 2 intersects B{p,£) in an open line segment. Let U = Cs\qi,q 2 - 
Since (T) does not contain a bnu, it follows that (M, oo) C U for some M < oo. Therefore, 
q' = qi + t{q 2 — qi) £ S' for a large enough integer t. Hence, there exists a point q £ S 
such that 'Ki{q) = q'■, so <I> is satisfiable. 

To prove the claim, let B = B(p,£) and let qi £ L\B. Consider the cone C = 
{qi + t{x — qi) \ X £ B,t > 0} and note that C contains B' := {g'l + d£~^{x — g^i) | x E 
B} = B{qi + 5£~^{p — qi),5). For a large enough positive constant 5, the set H' n L is 
non-empty. Let q 2 £ B' r\ L C. Then, the line through qi and q 2 intersects B in an 
open line segment. 

Finally, assume that (F) contains a bnu U. We may assume that U is not 0-valid: By 
Lemma uni (F) contains a non-empty unary relation T that is not 0-valid. Let c E Q be 
a non-zero constant such that U r\c-T ^ 0. IfC/flc-T contains more than one element, 
then it is a bnu that is not 0-valid. Otherwise, C/ H c ■ T is a constant unary relation, so 
(F) contains {!}. In this case, for a large enough constant c E Q, the relation U + c £ (F) 
is a bnu that is not 0 -valid. 

Let ri,r 2 E 1/ be two distinct points and let c E Q be a non-zero constant such that 
c-ri,C'r 2 E Z. Then, U' = c-U flZ is a bnu that excludes an interval and I7'n(—1,1) = 0. 
NP-hardness therefore follows from Lemma da □ 

7.2. The integer property. In this section, we will determine those semilinear constraint 
languages containing 7?+ for which knowing that there is a solution guarantees that there 
is an integer solution. We make the following definition. 

Definition 1. Let T he a constraint language over Q. We say that F has the integer 
property if every instance of CSP{T) has a solution if and only if it has an integer solution. 

The integer property can be used to infer tractability of certain semilinear constraint 
languages over Z. In particular, if F is a semilinear constraint language over Q that 
satisfies the integer property, then CSP(r) and CSP(F|z) are equivalent. To see that 
Fl^ is a semilinear constraint language over Z, take an arbitrary i? E F and let 99 be a 
quantifier-free definition of R over L/q[Z]. Then, p is also a quantifier-free definition of 
R\z over LIz[Z]. 

The following lemma shows that the integer property is preserved under pp-definitions. 
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Lemma 17. Let T be a constraint language over Q. IfV has the integer property, then 
so does (r). 

Proof. Let be an CSP-instance with relations Ri,..., Rk from (F), let be 

pp-definitions of Ri,..., Rk in F, and let be the CSF(r)-instance obtained from 'L 
by replacing each relation Ri by the quantifier-free part of tpi, and adding existential 
quantifiers for all auxiliary variables. If 'L has a rational solution, then 'L' has a rational 
solution, so 'h' has an integer solution. Note that the restriction of any solution of 'L' 
to Vars('I') is a solution to 'h. Therefore, the restriction of an integer solution of T' to 
Vars('I') is an integer solution to T, which proves the lemma. □ 

Let F denote a semilinear constraint language that contains R+. Observe that if {1} G 
(F), then CSF(r) cannot have the integer property since the following CSP(R+ U {!})- 
instance has the unique solution x = ^, y = 1: 

3x,y . X -h X = y A {l}y. 

Definition 2. Let T be a constraint language over Q. We say that F is scalable if the 
following holds: for each R gT and for each x = (xi,... ,Xk) G R, there exists a positive 
constant A such that {axi, . .., axk) G R, for all a > A. 

Clearly, scalable constraint languages cannot contain any unary constant relation {c} 
except when c = 0. Note that if F has endomorphisms e{x) = a ■ x for all rational 
a > A > 0, then F is indeed scalable. Inferring the existence of endomorphisms from the 
scalability property is, in general, not straightforward or even possible. The scalability 
property was originally defined slightly differently m but it is easy to verify that the two 
definitions coincide. 

The following result completely characterises the semilinear constraint languages that 
contain and have the integer property. 

Theorem 12. Let {ii+} C F C 5 Lq[Q] be a constraint language that is not 0-valid. Then, 
the following are equivalent: 

(1) T has the integer property. 

(2) every unary relation in (F) is either {0} or unbounded. 

(3) T does not satisfy (Poo)- 

(4) F is scalable. 

Proof. (HI) ^ ([2]). We show -i(l2]) ^ “■([I])- Suppose that Ti / {0} is a bounded unary 
relation in (F). By Lemma fOl there is a non-empty unary relation T 2 in (F) that is not 
0-valid. Therefore, for some c € Q, the unary relation U = Ti Ci c ■ T 2 in (F) is non-empty, 
bounded, and not 0-valid. Let k = 1+ [max(| sup U\, \ inf C/|)]. Consider the CSP instance 

3x, y . U{x) A k ■ y = X, 

and note that it has a solution: arbitrarily choose x G U and let y = x/k. However, it 
cannot have any integer solution since 0 and k was chosen such that k > |x|. Both U 
and the equation k ■ y = x are pp-definable in F, so the claim follows from Lemma fT71 

(ED ^ We show -dSD ^ -(ED- Assume that there exists a unary relation U in (F) 
containing a positive point and (M, 00 ) nU = 0 , for some M < 00 . If C/ is bounded, then 
“■ED follows immediately. Otherwise, by Lemma El there exists some M' < 00 such that 
{M', 00 ) nU = 0 and {— 00 ,—M') C U. By Lemma [H^2), there exists a bounded unary 
relation in ({i?_|_,17}) and, consequently, there exists such a relation in (F). 

ED ^ 0- We show -0 ^ -©. Arbitrarily choose an n-ary relation R G T such 
that R is not scalable. Arbitrarily choose a tuple p = (pi,... ,Pn) G R that witnesses that 
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R is not scalable, i.e., the set y = {y > 1 \ y ■ p € R} is unbounded. Consider the set 
U = {aGQ|a-p€ i?} and note that U is pp-definable in {i?, i?+} by Lemma [3j 

U{x) = 3yi, ...,yn.yi = x-piA---Ayn = x-pnA R{yi,.. .,yn)- 

Note that 1 £ U so U contains a positive point. Furthermore, since Y is unbounded, it 
follows from Lemma [2] that (M, oo) C Y for some M < oo. Hence, by definition of U, we 
have (M, oo) n [/ = 0, so F satisfies (Poo). 

(HD ^ (P). This implication is not difficult to deduce from the proof of Lemma 6 in m- 
We include an argument here for completeness. Assume that F is scalable and let $ be 
an arbitrary instance of CSP(r) with a solution x. Let i?i,... ,Rm be an enumeration 
of the atoms of <1> that contain a relation symbol from P. Since F is scalable, it follows 
that there exists a constant Ai such that ax satishes Ri{xi^,... ,Xi^) for all a > Ai. Let 
A = max{Ai,..., Am}- Then, ax satishes all atoms (including the equalities) of <1>, for all 
a> A. Therefore, if a is chosen to be a large enough common multiple of the denominators 
in x, then ax is an integral solution to <1?. □ 

As an immediate application of Theorem [T2] we give a complement to Theorem [TT] in 
the case when F violates (Pcxd)- 

Corollary 2. Let {i?+} C F C 5 Lq[Q] be a constraint language that violates (Poo)- The 
problem CSP{T U {Z}) is in P if 

• F is 0-valid; 

• (cone(r) U {{!}}) does not contain a bnu; or 

• cone(r) is essentially convex. 

Otherwise, CSP{T U Z) is NP-hard. 

Proof. If F is 0-valid, then F U {Z} is 0-valid, and hence in P. Otherwise, Theorem fT2] 
implies that F has the integer property. Therefore, every instance of CSP(r) has a solution 
if and only if it has an integer solution. It follows that CSP(r U {Z}) is polynomial-time 
equivalent to CSP(F). Since F violates (Poo), the result follows from Theorem [51 □ 

8. Discussion 

8.1. Generalisations. A natural goal, following the proof of Theorem |21 would be to 
determine the complexity of CSP(r) for an arbitrary semilinear constraint language T, 
i.e., when T does not necessarily contain R^. Below we indicate a few such attempts and 
the difficulties that accompany them. 

Consider Corollary [TJ Our main result. Theorem |21 generalises this by removing the 
assumption that {1} is in T. A natural question is then what happens if we instead 
remove the assumption that the addition relation needs to be in F. To this end, let SL^ 
denote the set of semilinear constraint languages such that {{!}} T F and {R+} ^ (F). A 
straightforward modification of the construction in Section 6.3 of Jonsson and Loow [16] 
gives the following: for every constraint language F' over a hnite domain, there exists a 
F G SL^ such that CSP(F') and CSP(F) are polynomial-time equivalent problems. Hence, 
a complete classihcation would give us a complete classihcation of finite-domain CSPs, 
and such a classihcation is a major open question within the CSP community UnidlKH]. 
We also observe that for every temporal constraint language (i.e., languages that are hrst- 
order dehnable in {<} over the rationals), there exists a F G SL^ such that CSP(F') 
and CSP(r) are polynomial-time equivalent problems. This follows from the fact that 
every temporal constraint language P' admits a polynomial-time reduction from CSP(r'U 
{{!}}) to CSP(F'): simply equate all variables appearing in {l}-constraints and note 
that any solution can be translated into a solution such that this variable is assigned the 
value 1. The complexity of temporal constraint languages is fully determined |1] and the 
polynomial-time solvable cases fall into nine different categories. The proof is complex 
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and it is based on the universal-algebraic approach for studying CSPs. We conclude that 
a complete classification of the languages in SL^ will require advanced techniques and will 
have to be conditioned on the classification of finite-domains CSPs. 

A smaller first step towards removing from Corollary [T] would be to only slightly 
relax the addition relation. Consider the affine addition relation A+ = {{a,b,c,d) G 
I a — b + c = d}. This relation can be viewed as a ‘relaxed’ variant of since A+ can 
be pp-defined in {R+} but not the other way round. Let P be a constraint language such 
that {A+,{c}} CPC 5 Lq[Q] for some c G Q. It is not hard to reduce the complexity 
classification for such constraint language to that of Theorem [2J 

Given a relation R C and a rational number c G Q, let i? -|- c denote the relation 
{(xi + c,..., Xk + c) I (xi,..., Xfc) G R}. For instance, A+ + c = A+. Similarly, we define 
r -|- c = {R -I- c I ii G r} for constraint languages P. Note that CSP(r) and CSP(r -|- c) 
are polynomial-time equivalent problems. 

Arbitrarily choose a constraint language {A+, {c}} CPC S'Lq[Q] and let P' = Pc). 
The problem CSP(P') is polynomial-time equivalent with CSP(P), A+ G P' and {0} G Pb 
The fact that A+ G P' and {0} G P' implies that R+ G (P') since R^{x,y,z) can be 
pp-defined by 

3w . {0}(r(;) A A+(x, w, y, z). 

Consequently, CSP(r') and CSP(P'U{i?+}) are polynomial-time equivalent problems. We 
conclude that CSP(r) is either in P or NP-complete by Theorem [21 

An interesting way forward would be to classify the complexity of CSP(r) for all {A_|_} C 
P C S'LqIQ]. Such a result would be a substantial generalisation of the results in Section 4 
of Bodirsky et al. [3]. Here, we see no obvious obstacles as in the case above for P G Sffi. 

8.2. The metaproblem. Theorem |2] shows that for every constraint language {i?+} C 
P C S'LqIQ], the problem CSP(P) is either in P or NP-complete. This makes the following 
computational problem (sometimes referred to as a metaproblem in the literature) relevant: 

Given a constraint language {R+} CPC S'LqIQ], is CSP(r) in P or NP-complete? 

We do not know the complexity of this problem and, in fact, it is not clear whether it is 
decidable or not. Interesting methods for tackling similar questions have been identified 
by, for instance, Bodirsky et al. [U] and Dumortier et al. BM- Bodirsky et al. analyse the 
decidability of abstract properties of constraint languages such as whether certain relations 
are pp-definable or not. Their results are based on a number of different techniques from 
model theory, universal algebra, Ramsey theory, and topological dynamics. Dumortier et 
al. [SI E] show that it is decidable whether a given first-order formula using the binary 
functions * and -|-, and the binary relation < over M with parameters from Q defines 
a semilinear relation. These results indicate that there are non-obvious properties of 
semilinear relations that may be relevant for proving (un)decidability of the metaproblem. 

Acknowledgements. The authors thank Manuel Bodirsky for suggesting the relation 
A+ as a relaxation of R+ (cf. Section [5TT|) . 
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